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Executive  Summary 

Research  Objective 

The  objective  of  this  research  is  to  develop  a  modeling  framework  for  simulating  the  insulative 
behavior  of  thermoplastic  polyurethane  elastomer  nanocomposites  (TPUNs)  for  solid  rocket 
motors  (SRMs)  and  the  creation  of  software  capable  of  predicting  the  properties  of  these  TPUNs. 

Research  Approach 

The  proposed  research  combines  numerical  modeling  and  experimental  characterization  of 
TPUNs  for  SRM  insulation.  The  TPUN  thermophyscial  properties  and  kinetic  parameters  are  to 
be  found  using  thermogravimetric  analysis,  Multi-Function  Sensor  (MFS)  gauges,  and  other 
thermal  analysis  techniques.  Numerical  techniques  for  small  scale  physics  from  molecular 
dynamics  (MD)  simulations  through  population  balance  methods  will  be  used  to  gain  insight  and 
intuition  on  the  interactions  between  the  polymeric  and  inorganic  materials.  These  insights  will 
help  us  to  develop  macroscopic  kinetic  models/descriptions  which  can  be  embedded  in 
macroscale  continuum  models.  These  models  will  be  first  exercised  against  the 
thermogravimetric  data  and  data  generated  from  the  MFS  gauges.  Subsequently,  small  scale 
experiments  using  the  LHMEL  facility  with  known  heating  conditions  and  well-controlled 
experiments  will  be  conducted  to  provide  us  with  data  to  further  validate  the  numerical  code.  The 
numerical  code  and  model  will  also  be  validated  with  mesoscale  experiments  using  AFRL  Pi-K 
and  Bates  SRM  firings.  Once  this  code  is  validated,  numerical  experiments  with  the  code  will 
help  us  to  design  optimized  compositions  (e.g.,  TPU  and/or  nanoparticles)  of  new  TPUNs. 

Scientific  Relevance 

The  mechanisms  through  which  the  addition  of  nanoparticles  modifies  the  decomposition  of 
thermoplastic  elastomers  are  not  fully  understood.  One  promising  method  for  understanding 
these  mechanisms  is  through  the  use  of  MD  simulations.  Because  of  the  time  scales  over  which 
MD  simulations  must  be  conducted,  it  has  not  been  obvious  how  to  incorporate  MD  findings  into 
continuum  type  simulations.  During  our  studies  it  is  anticipated  that  a  critical  scientific 
breakthrough  can  occur  if  we  are  able  to  develop  a  mechanistic  description  of  the  polymer- 
inorganic  material  coupling  process  that  can  reside  in  a  continuum  code.  One  approach  that  has 
merit  is  a  population  balance  (notional  particle)  description  of  the  polymer  and  inorganic 
material  distributions.  In  population  balance  treatment  of  aerosols,  kinetics  of  breakage  and 
coalescence  are  often  included  in  continuum  codes.  There  is  a  striking  analogy  between  polymer 
bond  breakage  and  subsequent  capture  by  nanoparticles  and  aerosol  aggregate  breakage  and 
subsequent  capture  by  a  scrubber  particle.  Our  ability  to  create  a  continuum  kinetic  description 
of  the  MD  physics  would  result  not  only  in  an  important  scientific  contribution  but  also  facilitate 
the  Air  Force’s  identification  and  design  of  new  and  future  compositions  for  SRM  insulation 
materials. 

Progress 

This  report  summarizes  the  first  four  months  of  our  research  effort  from  June  1  to  Sept.  30,  2006. 
For  this  period  of  performance,  we  summarized  a  simple  model  for  non-charring  ablation,  and 
then  used  that  model  to  perform  a  sensitivity  analysis  of  failure  time  based  on  four 
thermophysical  properties.  We  found  that  failure  time  is  most  sensitive  to  the  density  and  the 
heat  of  ablation  of  the  material.  It  was  also  found  that  failure  time  is  relatively  insensitive  to  the 
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thermal  conductivity  of  the  material.  These  results  indicate  that  improving  the  heat  of  ablation 
capacity  of  an  ablative  material  should  be  the  focus  of  future  research. 

We  also  reviewed  numerical  modeling  literature  for  thermal  protection  materials,  processed 
selective  TPUN  formulations  for  material  characterization,  and  calibrated  a  radiant  heat  panel 
apparatus  to  conduct  low  heat  flux  experiments  using  the  TPUNs.  A  user  proposal  entitled, 
“Thermophysical  Properties  Characterization  of  Thermoplastic  Elastomer  Nanocomposites”  was 
approved  by  Oak  Ridge  National  Laboratory/High  Temperature  Materials  Laboratory 
(ORNL/HTML).  This  allows  us  to  use  their  unique  material  characterization  equipment  and 
expertise  to  train  our  graduate  students  to  conduct  thermophysical  property  measurements  of  our 
TPUNs. 

Personnel  Supported 

Morgan  Bruns  PhD  student,  Mechanical  Engineering 

Dave  Ho  MS  student,  Mechanical  Engineering 

Ofodike  A.  Ezekoye,  Ph.D,  P.E.  Professor,  Mechanical  Engineering 

Joseph  H.  Koo,  Sc.D.  Senior  Research  Fellow,  Mechanical  Engineering 
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1.  Introduction 

Insulation  materials  are  required  to  protect  solid  rocket  motor  (SRM)  cases  against  hostile 
combustion  gaseous  products  generated  in  the  rocket  chamber.  Various  internal  system 
components  must  be  protected  from  extreme  flow  temperatures  of  1,000  to  4,000°C  and  highly 
abrasive  particles  ejected  at  velocities  greater  than  lkm/s.  Kevlar®- filled  ethylene-propylene- 
diene  rubber  (EPDM)  is  the  baseline  insulation  material  for  solid  rocket  motor  cases.  A  novel 
class  of  insulation  materials  was  developed  by  the  Air  Force  Research  Laboratory  that  is  lighter, 
exhibits  better  erosion  and  insulation  characteristics  and  possesses  a  more  cost-effective 
manufacturing  process  than  the  current  baseline  material.  ‘‘5  We  propose  to  develop  a 
comprehensive  numerical  and  experimental  research  program  to  understand  the  ablation  and  heat 
transfer  characteristics  of  this  new  class  of  thermoplastic  elastomer  nanocomposites  to  facilitate 
the  Air  Force’s  identification  and  design  efforts  for  new  and  future  components  of  SRM 
insulation  materials. 

Koo  et  al.6~i  have  investigated  different  types  of  nanoparticles  besides  MMT  nanoclays  with 
phenolic  resole  resin  for  nozzle  ablative  material.  These  nanostructured  materials  exploit  the 
ablation  resistance  of  both  resin  and  nanoparticles.  The  objective  of  this  study  was  the 
development  of  a  new  class  of  nanostructured  materials  that  are  lighter  and  have  better  erosion 
and  insulation  characteristics  than  current  ablatives  for  rocket  nozzle  assembly.  Borden 
Chemicals’  SC- 1008,  a  phenolic  resole,  was  selected  as  the  resin.  Several  nanoparticles  such  as 
Southern  Clay  Products  (SCP)  nanoclay,  Applied  Sciences  carbon  nanofibers  (CNFs),  and 
Hybrid  Plastics  POSS®  were  evaluated  with  the  SC- 1008  phenolic  resin.  Based  on  wide-angle 
x-ray  diffraction  (WAXD)  and  transmission  electron  microscopy  (TEM),  selected  nanoparticles 
were  identified  for  further  study.  Cytec  Engineered  Materials  (CEM)  fabricated  several  MX- 
4926  alternates  by  replacing  the  carbon  black  filler  with  the  selected  nanoparticles  for  ablation 
testing.  These  polymer  nanocomposites  demonstrated  significantly  improved  ablation  and 
insulative  performance  versus  the  baseline  ablative.’1"* 

A  subscale  solid  rocket  Char  motor  developed  by  AFRL/PRSM  was  used  to  study  the  ablation 
characteristics  of  the  insulation  materials.1  Testing  insulative  or  ablative  materials  in  full-scale 
firings  is  not  only  expensive,  but  also  requires  lengthy  planning  and  limited  exposure  time.  Thus 
reducing  the  number  of  samples  and  the  variety  of  test  conditions  available  is  a  worthwhile 
endeavor.  The  Char  motor  has  been  shown  to  be  a  cost-effective,  “fast  turnaround”  device  to 
evaluate  different  insulative  materials  under  identical  conditions  for  initial  material  screening 
and  development.  Physical,  mechanical,  and  thermophysical  properties  of  these  novel  polymer 
nanocomposites  have  been  determined  and  compared  with  baseline  materials." 

2.  Technical  Program  Summary 

During  the  first  year  of  the  three-year  program,  processing  of  selective  TPUNs  in  collaboration 
with  AFRL/PRSM  will  be  conducted.  Multi-Function  Sensor  (MFS)  gauges  using  TPUNs  will 
be  designed  and  fabricated.  Characterization  experiments  to  understand  the  degradation  physics 
of  these  TPUNs  will  be  conducted.  TPUNs  thermophysical  properties  and  kinetic  parameter  will 
be  determined.  Molecular  dynamics  (MD)  simulation  of  representative  systems  will  be 
evaluated  to  gain  insight  into  the  role  of  nanoparticles  on  polymer  chain  stability. 
Phenomenological  models  using  experimental  data,  MD  results,  and  literature  will  be  developed. 
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Requirements  of  TPUN  thermal  degradation  computational  tools  will  be  determined.  Research 
progress  will  be  documented  in  a  first  year  progress  report. 

Under  year-two  funding,  a  numerical  code  will  be  developed.  Validation  experiments  with  small 
scale  degradation  physics  using  the  AFRL  LHMEL  facility  will  be  performed.  LHMEL  data 
will  be  used  to  test  numerical  code.  Validation  experiments  at  mesoscales  will  be  carried  out 
using  AFRL/PRSM’s  2-lb  Pi-K  SRM  firings  in  collaboration  with  AFRL  personnel.  The 
emerging  new  code  will  be  used  to  conduct  sensitivity  tests.  We  will  participate  in  designing  the 
AFRL  15-lb  Bates  motor  case  with  heat  flux  and  MFS  gauges  to  collect  ablation  data  for  the  new 
TPUNs.  Processing  of  new  TPUNs  in  collaboration  with  AFRL  will  be  accomplished. 
Thermophyscial  properties  and  kinetic  parameters  for  the  new  TPUNs  will  be  measured.  MFS 
gauge  for  new  TPUNs  will  be  fabricated.  Research  progress  will  be  documented  in  a  second 
year  progress  report. 

Year-three  funding  will  be  used  to  refine  the  numerical  code  using  validation  experiments  of 
small  scale  degradation  physics  of  the  new  TPUNs  using  the  LHMEL  facility.  Validation 
experiments  of  the  mesoscale  15-lb  Bates  SRM  with  collaboration  with  AFRL  will  be  conducted 
and  data  analyzed.  A  user’s  manual  for  the  new  code  will  be  developed  with  sample  cases. 
Research  progress  for  the  entire  program  will  be  documented  in  a  final  report. 

3.0  Program  Work  Plan 

A  chronological  list  of  program  tasks  by  program  years  and  a  program  schedule  is  listed  in  Table 
1.  Key  research  challenges  involved  in  the  proposed  research  are  discussed  below. 
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Table  1  Milestone  for  the  Experimental  and  Numerical  Characterization  of  Polymer 
Nanocomposites  for  Solid  Rocket  Motor  Internal  Insulation  Research  Program 


Task  Description 

First  Year  Tasks 

Fisrt  Year 

Q1  Q2  Q3  Q4 

Second  Year 

Q5  Q6  Q7  Q8 

Third  Year 

Q9  Q10Q11  Q12 

Task  1  -  Process  Selective  TPUNs 

Task  2  -  Design  &  Fabricate  MFS  Gauge  with  TPUN 
Task  3  -  Conduct  Characterization  Expts 

Task  4  -  Determine  TPUN  Properties  &  Kinetics 

Task  5  -  Develop  Phenomenological  Models 

Task  6  -  Determine  Reqs  of  TPUNTD  Computational  T 
Task  7  -  First  Year  Progress  Report/Technical  Papers 

X 

X  X 

XXX 

X  X 

X  X  X  X 

ool  X  X 

X 

Second  Year  Tasks 

Task  8  -  Implement  Numerical  Code 

Task  9  -  Perform  Validation  Expts  Using  LHMEL 

Task  10  -  Test  Code  with  LHMEL  Data 

Task  11  -  Perform  Validation  Expts  Using  Pi-K  Motor 
Task  12  -  Conduct  Sensitivity  Test  Using  New  Code 
Task  13  -  Design  15-lb  Bates  Motor  with  Different  Gau< 
Task  14  -  Process  New  TPUNs 

Task  15  -  Determine  New  TPUN  Properties  and  Kinetic 
Task  16  -  Fabricate  MFS  Gauge  for  New  TPUNs 

Task  17  -  Second  Year  Progress  Report/Technical  Pap 

]es 

1 

:s 

1 

>ers 

X  X 

X 

X  X 

X  X 

X  X 

X  X  X  X 

X  X 

X  X 

X  X 

X 

Third  Year  Tasks 

Task  18  -  Perform  Validation  Expts  Using  LHMEL 

Task  19  -  Perform  Validation  Expts  Using  Bates  Motor 
Task  20  -  Generate  User's  Manual 

Task  21  -  Final  Program  Report/Technical  Papers 

X  X 

X  X 

X  X 

X 

3.1  First  Year  Program  Tasks 

■  Task  1  Process  Selective  TPUNs  in  Collaboration  with  AFRL/PRSM,  Edwards  AFB  - 
Based  on  previous  studies1'  PELLETHANE™  2102-90A  with  5%  Cloisite  30B  and 
PELLETHANE™  2102-90A  with  15%  PR-19-PS  CNF  TPUNs  have  the  best  ablation 
resistant  and  mechanical  properties.  These  two  materials  are  most  likely  our  prime 
candidates  for  our  initial  candidates.  The  Desmopan®  DP  6065A-Cloisite  30B  and  the 
DP  6065A-PR-19-PS  CNF  are  also  potential  candidates.  These  PNMs  will  be  provided 
to  our  program  by  AFRL/PRSM. 

■  Task  2  Design  and  Fabricate  Multi-Function  Sensor  (MFS)  Gauge  Using  TPUNs  - 
Several  types  of  MFS  gauges  will  be  designed  with  NANMAC  engineers  and  fabricated 
for  UT  materials  characterization  and  LHMEL  testing  in  Yr.  2. 

■  Task  3  Conduct  Characterization  Experiments  (Degradation  Physics)  @UT  - 
Characterization  experiments  will  be  conducted  in  the  UT  laboratory  to  understand  the 
ablation  and  heat  transfer  characteristics  of  these  TPUNs.  We  will  more  clearly  explore 
the  physics  of  thermoplastic  polyurethane  elastomer  nanocomposite  thermal  degradation 
(TPUNTD).  In  particular  we  will  conduct  small  scale,  well  controlled,  high  heat  flux 
thermal  degradation  experiments.  Using  observations  (visual  and  post  processed 
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microscopy)  and  data  as  generated  by  in-situ  thermocouples  we  will  develop  simple 
phenomenological  models  for  the  degradation  process. 

■  Task  4  Determine  TPUN  Properties  and  Kinetic  Parameters/Detailed  Measurements  of 
thermophysical  properties  -  Based  on  degradation  physics  a  list  of  thermophysical 
properties  based  on  Task  3  will  be  conducted  at  UT  and  Oak  Ridge  National 
Laboratory/High  Temperature  Materials  Laboratory  by  UT  graduate  students. 

■  Task  5  Develop  Phenomenological  Models  using  Experimental  Data,  Molecular 
Dynamics  Simulations,  &  Literature  -  phenomenological  models  will  be  developed  to 
simulate  material  behavior  under  SRM  chamber  environment.  We  will  explore  the 
literature  for  potential  approaches  to  develop  these  models.  The  phenomenological 
model  will  be  the  basis  for  development  of  our  TPUNTD  computational  tool. 

■  Task  6  Determine  Requirements  of  TPUNTD  Computational  Tools  -  The  functional 
requirements  for  the  tool  will  be  determined  based  upon  the  level  of  sophistication  that 
can  be  achieved  by  the  phenomenological  models.  As  an  example,  the  tools  can  be  used 
in  one  of  three  ways.  In  order  of  simplest  uses  to  more  sophisticated  use  we  imagine  that 
the  tools  can  be  used  for:  scaling  up  of  reduced  scale  tests,  prediction  of  performance  of 
new  TPUN  materials,  and  design  and  optimization  of  material  parameters  for  thermal 
function. 

■  Task  7  Document  activities  and  results  into  a  first  year  progress  report. 

3.2  Second  Year  Program  Tasks 

■  Task  8  Implement  Numerical  Code  -  Once  the  phenomenological  model  has  been 
constructed,  it  will  be  necessary  to  determine  the  proper  numerical  analysis  methods  to 
apply  to  the  problem.  It  is  anticipated  that  we  will  initially  develop  a  one-dimensional 
approximation.  We  will  likely  use  a  high  order  finite  difference  method  for  the 
discretization.  From  the  phenomenological  model  TPUNTD  we  will  know  which 
property  data  will  be  required.  As  such,  we  will  perform  appropriate  experimental 
characterization  of  the  materials  to  generate  this  property  data. 

■  Task  9  Perform  Validation  Expts  (Small  Scale  Degradation  Physics)  Using  LHMEL 
(Round  1)  -  The  LHMEL  facility  will  be  used  to  collect  data  to  test  our  numerical  model. 
LHMEL  I  will  be  selected  for  our  study.  Test  conditions  (laser  power  and  test  duration) 
will  be  determined  later  with  guidance  from  AFRL/PRSB,  Edwards  AFB  personnel. 

■  Task  10  Test  code  with  LHMEL  Data  -  LHMEL  test  data  will  be  used  to  verify  the  code. 

■  Task  11  Perform  Validation  Expts  (Meso  Scale  AFRL/Edwards  2-lb  Pi-K  Solid  Rocket 
Motor)  in  Collaborate  with  AFRL/PRSM,  Edwards  AFB.  This  task  will  enable  us  to  test 
TPUNs  in  a  solid  rocket  motor  chamber  environment  and  collect  data  to  verify  the 
numerical  model. 

■  Task  12  Conduct  Sensitivity  Tests  using  New  Code. 

■  Task  13  Design  15-lb  Bates  SRM  Case  in  Collaborate  with  AFRL/PRSM  using  Heat 
Flux  and  MFS  Gauges  to  Collect  Ablation  Data  with  New  TPUNs  (Round  1).  Heat  flux 
and  MFS  gauges  with  New  TPUNs  will  be  installed  in  the  Bates  motor  to  measure  the 
heat  flux  at  different  locations  of  the  cone  section  in  the  motor  chamber.  Ablation  and 
heat  transfer  data  will  be  collected  to  test  our  new  code. 
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■  Task  14  Process  New  TPUNs  in  Collaborate  with  AFRL/PRSM  -  New  types  of  TPU 
with  different  nanoparticles  will  be  investigated. 

■  Task  15  Determine  Properties  and  Kinetic  Parameters  for  New  TPUNs  -  Thermophysical 
properties  of  the  new  TPUNs  will  be  characterized  at  Oak  Ridge  National 
Laboratory/High  Temperature  Materials  Laboratory  for  model  verification. 

■  Task  16  Fabricate  Multi-Function  Sensor  (MFS)  Gauge  for  New  TPUNs  -  MFS  gauges 
will  be  fabricated  for  LHMEL  and  SRM  testing. 

■  Task  17  Document  activities  and  results  into  a  second  year  progress  report. 

3.3  Third  Year  Program  Tasks 

■  Task  18  Perform  Validation  Expts  (Small  Scale  Degradation  Physics)  Using  LHMEL 
(Round  2)  -  Collect  ablation  and  heat  transfer  data  for  code  validation. 

■  Task  19  Perform  Validation  Expts  [Meso  Scale  AFRL/Edwards  15-lb  Bates  SRM  (Round 
2)]  in  Collaborate  with  AFRL/PRSM,  Edwards  AFB  -  Collect  ablation  and  heat  transfer 
data  for  code  validation 

■  Task  20  Generate  User’s  Manual  -  a  user  friendly  manual  will  be  developed  with  sample 
cases. 

■  Task  21  Document  activities  and  results  into  final  program  report. 

4.0  Description  of  Polymer  Nanocomposites 

4.1  Thermoplastic  Elastomer 

Thermoplastic  elastomer  (TPE)  is  the  identifying  term  for  a  family  of  rubberlike  materials  that, 
unlike  conventional  vulcanized  rubbers,  can  be  processed  and  recycled  like  thermoplastic 
materials.  One  example  of  a  thermoplastic  elastomer  is  PELLETHANE™  2102-90A. 
PELLETHANE™  2102-90A  is  a  thermoplastic  polyurethane  elastomer  (TPU)  manufactured  by 
Dow  Chemicalf  Its  typical  applications  include  seals,  gaskets,  belting,  animal  ID  tags,  and  other 
fabricated  products.  It  not  only  has  excellent  resistance  to  fuels  and  oils,  but  also  has  good 
hydrolytic  stability. 

Another  example  of  TPU  is  Desmopan®  DP  6065A,  a  special  ether  TPU  manufactured  by  Bayer 
Material  Science.1 '  Typical  uses  are  automotive  instrument  panels,  caster  wheels,  power  tools, 
sporting  goods,  medical  devices,  and  a  variety  of  extruded  film,  sheet  and  profile  applications. 
This  TPU  is  the  new  candidate  that  AFRL/PRSM  has  selected  to  be  modified  with  different 
nanostructures. 

4.2  Montmorillonite  Nanoclays 

Exfoliation  of  surface  treated  montmorillonite  (MMT)  clays  into  a  polymeric  material  have  been 
shown  to  improve  mechanical  properties,  barrier  performance,  and  application  processing.11, 12 
Achieving  exfoliation  of  organomontmorillonite  in  various  continuous  phases  is  a  function  of  the 
surface  treatment  of  the  MMT  clays  and  the  mixing  efficiency  of  the  dispersing  protocol. 
Surface  treatment  of  MMT  is  classically  accomplished  with  the  exchange  of  inorganic 
counterions,  e.g.,  sodium,  etc.,  with  quaternary  ammonium  ions.  The  choice  of  the  quaternary 
ammonium  ion  is  a  function  of  the  hydrophilic/hydrophobic  nature  of  the  continuous  phase,  in 
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this  case  a  thermoplastic  polyurethane  elastomer  polymer.  The  quaternary  ammonium  ion  acts 
as  a  compatibilizing  agent  between  the  interface  of  MMT  and  the  continuous  phase  of  the  TPU. 

The  Cloisite®  30B  is  a  surface  treated  montmorillonite  (Tallow  bishydroxyethyl  methyl, 
T(EOH)2M)  manufactured  by  SCP.1'  It  is  an  additive  for  plastics  to  improve  various  plastic 
physical  properties,  such  as  modulus,  heat  distortion  temperature,  and  permeability.  It  has  a 
ternary  ammonium  salt  modifier  at  90  meq/lOOg  clay  and  is  off-white  in  color. 

In  an  earlier  study  by  Koo  et  al.  ‘  Cloisite®  30B  in  loadings  of  5,  10,  and  15  wt%  was  dispersed 
in  a  phenolic  resole  (SC- 1008).  The  TEM  analyses  allowed  us  to  determine  the  degree  of 
dispersion/exfoliation  of  the  nanoclay  before  committing  to  a  20-lb  run  at  CEM  of  these 
nanoparticle/resin  mixtures  to  make  prepregs.  The  resulting  prepregs  were  transformed  into 
molded  components  for  solid  rocket  nozzle  assembly.  This  proved  to  be  a  very  cost-effective 
and  efficient  technique  for  screening  different  formulations.  This  nanoclay  is  quite  versatile  and 
has  been  effective  for  phenolic,6"’'  cyanate  ester,  4  ‘  and  epoxy  resin  systems.1, 

4,3  Carbon  Nanofibers  (CNFs) 

CNFs  are  a  form  of  vapor-grown  carbon  fiber,  which  is  a  discontinuous  graphitic  filament 
produced  in  the  gas  phase  from  the  pyrolysis  of  hydrocarbons.1  In  terms  of  physical  size, 
performance  improvement,  and  product  cost,  CNFs  complete  a  continuum  bounded  by  carbon 
black,  fullerenes,  and  single-wall  to  multi-wall  carbon  nanotubes  on  one  end  and  continuous 
carbon  fiber  on  the  other.  CNFs  are  able  to  combine  many  of  the  advantages  of  these  other 
forms  of  carbon  for  reinforcement  in  engineered  polymers.  The  CNFs  have  transport  and 
mechanical  properties  that  approach  the  theoretical  values  of  single  crystal  graphite,  similar  to 
the  fullerenes,  but  can  be  made  in  high  volumes  at  low  cost  —  ultimately  lower  than  conventional 
carbon  fibers.  In  equivalent  production  volumes,  CNFs  are  projected  to  have  a  cost  comparable 
to  E-glass  on  a  per-pound  basis,  yet  possess  properties  which  far  exceed  those  of  glass  and  are 
equal  to,  or  exceed  those  of  much  more  costly  commercial  carbon  fiber.  Maruyama  and  Alam 
published  an  excellent  review  of  carbon  nanotubes  and  nanofibers  in  composite  materials."1 

The  CNFs  are  fabricated  by  Applied  Sciences  Inc./Pyrograf®  Products  by  pyrolytic 
decomposition  of  methane  in  the  presence  of  iron-based  catalyst  particles  at  temperatures  above 
900°C.  Pyrograf®-III  is  a  patented,  very  fine,  highly  graphitic,  CNFs.  Pyrograf®-III  is  available 
in  diameters  ranging  from  50  to  200  nm  and  a  length  of  50  to  100  pm.  Therefore,  CNFs  are 
much  smaller  than  conventional  continuous  or  milled  carbon  fibers  (5  to  10  pm)  but  significantly 
larger  than  carbon  nanotubes  (1  to  10  nm).  Compared  to  PAN  and  pitch-based  carbon  fiber,  the 
morphology  of  CNFs  is  unique  in  that  there  exists  far  less  impurities  in  the  filament,  providing 
for  graphitic  and  turbostatic  graphite  structures.  In  addition,  the  graphene  planes  are  more 
preferentially  oriented  around  the  axis  of  the  fiber.  Consequences  of  the  circumferential 
orientation  of  high  purity  graphene  planes  are  a  lack  of  cross-linking  between  the  graphene 
layers,  and  a  relative  lack  of  active  sites  on  the  fiber  surface,  making  it  more  resistant  to 
oxidation,  and  less  reactive  for  bonding  to  matrix  materials.  Also  in  contrast  to  carbon  fiber 
derived  from  PAN  or  pitch  precursors,  CNFs  are  produced  only  in  a  discontinuous  form,  where 
the  length  of  the  fiber  can  be  varied  from  about  100  pm  to  several  cm,  and  the  diameter  is  of  the 
order  of  100  nm. 

The  CNFs  have  exceptional  mechanical  and  transport  properties,  providing  them  with  excellent 
potential  as  an  engineering  material.  Such  fibers  consist  almost  exclusively  of  CVD  carbon, 
which  is  less  graphitic  and  more  defective  than  the  catalytically  grown  carbon  core  that 
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constitutes  the  CNFs.  Thus,  the  properties  listed  in  the  table  represent  an  estimate  for  the 
properties  of  nanofibers.  The  two  CNFs  used  in  the  study  are  PR-19-PS  and  PR-24-PS. 

4.4  Ablation  of  TPUNs 

Eight  of  the  thirteen  TPU  blends  were  tested  with  the  baseline  Kevlar®-filled  EPDM  using  the 
Air  Force  Research  Laboratory  Char  motor  using  aluminized  propellants.  Table  2  shows  a 
summary  of  TPU  blends,  wt%  of  filler,  theoretical  and  actual  density,  and  ablation  ranking.1’2 
Ranking  of  1  indicates  the  best  ablation  resistant  material.  The  5%  Cloisite®  30B  TPUN  and 
15%  PR-19-PS  CNF  TPUN  are  ranked  highest  for  their  respective  family  of  nanoparticles. 

Table  2  Summary  of  Densities  and  Ablation  Ranking  of  Kevlar®-EPDM  and  TPUNs 


Polymer 

Matrix 

Filler 

Wt% 

Theoretical 

Density 

(g/cc) 

Actual 

Density 

(g/cc) 

Ablation 

Ranking 

Kevlar®- 

Kevlar® 

12 

0.95 

0.97 

8 

EPDM 

2 102-90 A 

Cloisite®  30B 

5 

1.20 

1.21 

1 

2102-90A 

Cloisite®  30B 

7.5 

1.21 

1.22 

3 

2 102-90 A 

Cloisite®  30B 

10 

1.23 

1.23 

2 

2 102-90 A 

PR-24-PS  CNF 

5 

1.20 

1.18 

10 

2 102-90 A 

PR-24-PS  CNF 

15 

1.22 

1.26 

6 

2 102-90 A 

PR-24-PS  CNF 

25 

1.25 

1.30 

9 

2102-90A 

PR-19-PS  CNF 

5 

1.20 

1.21 

4 

2102-90A 

PR-19-PS  CNF 

15 

1.25 

1.27 

5 

2 102-90 A 

PR-19-PS  CNF 

25 

1.22 

1.26 

6 

4.5  Microsctructural  Analyses  of  TPUNs 

Post-test  specimens  of  Kevlar®-filled  EPDM,  15%  PR-24-PS  CNFs:  85%  2102-90A  TPU,  and 
7.5%  Cloisite®  30B:  92.5%  2102-90A  TPU  were  analyzed  using  SEM.  The  characterization  of 
the  nanostructured  materials  using  SEM  enables  us  to  understand  the  behavior  of  these  materials. 
SEM  images  were  produced  to  examine  the  surface  and  x-sectional  views  of  these 
nanostructured  materials. 

Figures  1(a)  to  1(c)  show  the  SEM  images  of  the  surface  views  of  Kevlar®-filled  EPDM,  15% 
PR-24-PS  CNFs:  85%  2102-90A  TPU,  and  7.5%  Cloisite®  30B:  92.5%  2102-90A  TPU 
specimens,  respectively.  The  scale  bars  are  10  pm  in  the  SEM  micrographs.  The  Kevlar®-filled 
EPDM  material  (Figure  1(a) )  was  exposed  to  non-aluminized  propellant  combustion  gaseous 
products.  Complementary  15%  PR-24-PS  CNFs:  85%  2102-90A  TPU  post-test  specimen  from 
the  same  motor  firing  is  shown  in  Figure  1  (b).  Figure  1(c)  shows  the  SEM  image  of  7.5% 
Cloisite®  30B:  92.5%  2102-90A  TPU  post-test  specimens  from  a  different  motor  firing.  These 
materials  were  exposed  to  aluminized  propellant  combustion  gaseous  products.  It  was  noted  that 
the  surface  behavior  of  the  charred  Kevlar®-EPDM,  nanoclay  TPUN,  and  CNF  TPUN  are  very 
different  as  shown  in  Figure  1.  Upon  analysis  of  the  control  sample,  it  was  observed  that  the 
EPDM  rubber  was  totally  burned  forming  a  small  amount  of  char,  which  was  easily  removed 
from  the  remaining  EPDM  rubber.  In  contrast,  the  testing  of  the  TPUNs  resulted  in  a  larger 
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amount  of  char  that  was  difficult  to  remove  from  the  remaining  TPU.  The  lower  ablation  rate 


coupled  with  the  more  structurally  stable  char  indicates  higher  performance  insulation. 


(a)  (b)  (c) 

Figures  1  (a)  to  1(c)  (a)  SEM  micrographs  of  surface  views  of  Kevlar-EPDM,  (b)  15%  PR-24- 
PS  CNFs:  85%  2102-90A  TPU,  and  (c)  7.5%  Cloisite  30B:  92.5%  2102-90A  TPU;  scale  bars  are 

10  pm. 

5.0  Discussion  of  Results 

This  annual  report  summarizes  the  first  four  months  of  our  research  progress  during  the  period  of 
June  1  to  Sept.  30,  2006.  This  section  will  be  divided  into  four  sections:  (a)  Review  paper,  (b) 
Experimental  progress,  (c)  Numerical  progress,  and  (d)  Publications. 

5.1  Review  Paper 

Thermal  protection  materials  are  required  to  protect  structural  components  of  space  vehicles 
during  the  re-entry  stage,  missile  launching  systems,  and  solid  rocket  motors.  A  thorough 
literature  survey  was  conducted  to  review  the  numerical  and  experimental  characterization  of 
these  thermal  protection  materials  for  different  military  and  aerospace  applications.  The 
literature  survey  is  grouped  into:  (a)  numerical  modeling,  (b)  experimental  testing,  (c)  materials 
thermophysical  properties  characterization,  and  (d)  advanced  nozzle  throat  material.  The 
numerical  modeling  portion  was  completed  and  will  be  summarized  in  this  section,  the  other 
portions  of  the  literature  survey  are  still  in  progress. 

In  this  section,  only  the  conclusions  of  the  numerical  modeling  literature  will  be  discussed 
detailed  discussion  can  be  found  in  Ref.  22.  The  experimental  testing,  materials  thermophysical 
property  characterization,  and  advanced  nozzle  throat  material  will  be  presented  in  subsequent 
reviews.  Numerical  modeling  literature  is  subdivided  into  four  groups  based  on  material 
modeling  and  the  computational  fluid  dynamic  (CFD)  analyses  that  were  used  for  different 
applications.  The  different  application  groups  are  described  such  as:  (a)  thermal  protection 
systems  (TPS)  for  spacecrafts,  (b)  ablative  material  for  missile  launching  systems,  (c)  internal 
insulation  for  solid  rocket  motors,  and  (d)  nozzle  assembly  for  solid  rocket  motors. 

Reviews  of  numerical  models  and  CFD  analyses  are  presented  in  chronological  order  in  Ref.  22. 
The  following  conclusions  are  presented: 

1.  Computational  fluid  dynamic  analyses  are  usually  coupled  with  material  response 
modeling  to  provide  a  realistic  solution  for  each  application.  This  coupling  is  well 
developed  for  spacecraft  and  missile  areas  in  one-  and  two-dimensional  analyses.  The 
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Charring  Material  Thermal  Response  and  Ablation  Program  (CMA)^  is  the  industry 
standard  code  for  this  type  of  application.  The  development  and  validation  of  the  three- 
dimensional  CFD  and  CMA  analyses  are  in  progress. 

2.  Limited  data  can  be  obtained  in  the  open  literature  as  it  relates  to  coupling  CFD  analyses 
with  material  modeling  for  the  missile  launching  system  applications. 

3.  There  are  a  few  studies  in  coupling  CFD  analyses  with  material  modeling  in  the  solid 
rocket  motor  insulation  application.  The  CFD  analyses  of  the  SRM  flowfield  is  more 
advanced  than  the  material  modeling  of  insulation  material  such  as  Kevlar®  filled 
EPDM.  Treatment  of  SRM  insulation  modeling  is  mostly  empirical.  Material  modeling 
effort  is  warranted. 

4.  There  are  current  activities  in  the  technical  community  in  conducting  advanced  CFD 
analyses  coupling  with  nozzle  material  response  modeling. 

5.  Material  properties  characterization  is  urgently  needed  for  new  materials  to  facilitate  their 
use  in  the  industry  standard  codes.  A  review  of  the  material  properties  characterization 
literature  is  in  progress  by  the  authors.24 

5.2  Experimental  Progress 

TPUN  Processing 

Four  loading  levels  in  2.5,  5,  7.5,  and  10  wt%  of  Cloisite  30B  were  incorporated  into  the  2102- 
90A  TPU  to  create  polymer  nanocomposites  using  twin  screw  extrusion  at  2Ll  Century 
Polymers.  Similarly  three  levels  in  10,  15,  and  20  wt%  of  PR-19-PS  CNFs  were  also 
incorporated  into  the  2102-90A  TPU  using  twin  screw  extrusion  processing.  This  task  was  done 
in  conjunction  with  an  Aerojet/AFRL  program.  TEM  analyses  will  be  conducted  for  all 
formulations  to  ensure  good  dispersion  of  these  nanoparticles  in  the  polymer.  Selective 
candidates  from  this  material  matrix  will  be  tested  using  the  Simulated  Solid  Rocket  Motor 
(SSRM)  at  Texas  Sate  University-San  Marcos,  TX  at  high  heat  flux  such  as  14  MW/m2.  The  test 
series  are  scheduled  to  be  performed  in  Oct. 

TPUN  Characterization 

TGA  and  DSC  experiments  at  heating  rates  of  5,  10,  and  20°C/min  are  in  progress  to  study  the 
thermal  properties  of  these  TPUNs  at  the  UT-Austin.  Besides  thermal  characterization, 
microstructural  analyses  of  these  pre-  and  post-test  TPUN  specimens  will  be  conducted  using 
SEM. 

A  research  proposal  entitled,  “Thermophysical  Properties  Characterization  of  Thermoplastic 
Elastomer  Nanocomposites’’  was  submitted  to  Oak  Ridge  National  Laboratory/High 
Temperature  Materials  Laboratory.  Under  DOE  sponsorship  ORNL/HTML  has  unique 
thermophysical  properties  characterization  equipment  and  staff  to  train  our  graduate  students  to 
conduct  thermophysical  property  measurements.  Our  proposal  was  approved  by  ORNL/HTML 
recently.  One  of  our  graduate  students  and  J.  Koo  will  be  visiting  ORNL/HTML  to  conduct  the 
following  thermophyscial  property  measurements: 

•  Temperature-dependent  measurements  of  specific  heat  for  virgin  and  charred  TPUN-clay 
and  TPUN-CNFs 
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•  Temperature-dependent  measurements  of  thermal  diffusivity  for  virgin  and  charred 
TPUN-clay  and  TPUN-CNFs 

•  Temperature-dependent  measurements  of  thermal  expansion  for  virgin  and  charred 
TPUN-clay  and  TPUN-CNFs 

Most  likely  the  5  wt%  Cloisite  30B  TPUN  and  the  15  wt%  PR-19-PS  CNFs  TPU  will  be  selected 
for  the  detailed  thermophysical  properties  characterization.  These  temperature-dependent 
material  properties  will  be  used  for  the  development  of  numerical  models  for  this  program. 

Radiant  Heat  Panel 

A  quartz  radiant  heat  panel  that  can  achieve  heat  flux  of  82  kW/m  was  set  up  and  calibrated  to 
conduct  low  heat  flux  fire  simulations  in  our  Combustion  Laboratory.  Figures  2  (a)  and  (b) 
show  the  heat  flux  gauge  used  to  calibrate  the  heat  flux  of  the  radiant  panel  and  radiant  panel 
experimental  setup,  respectively. 


(a)  (b) 

Figures  2  (a)  and  (b)  (a)  Heat  flux  gauge  mounted  on  carbon/carbon  composite  to  calibrate  heat 
flux  of  the  quartz  radiant  panel  and  (b)  Radiant  heating  panel  with  sample  exposing  to  the  red  hot 
radiant  panel. 

Figure  3  shows  the  calibration  curve  of  the  radiant  panel  versus  distance  of  the  test  specimen 
away  from  the  radiant  panel.  Several  TPUNs  have  been  tested  to  check  out  the  system.  This 
experimental  setup  will  enable  us  to  characterize  the  TPUNs  at  known  heat  fluxes  and  SEM 
microstructural  analyses  will  enable  us  to  learn  about  material  behavior  under  simulated  fire 
conditions. 
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Figure  3  Calibration  of  heat  flux  versus  distance  from  radiant  heat  panel. 

5.3  Numerical  Progress 

Progress  has  been  made  on  three  fronts  in  the  numerical  and  modeling  portion  of  this  project. 
We  have  begun  developing  a  simulation  method  for  the  intermediate  physics  part  of  the  project 
where  the  population  balance  modeling  fits  most  naturally.  We  have  been  working  with  the 
traditional  continuum  based  models  (Hero-2D  and  ACE'  )  to  better  understand  the  limitations 
of  those  models.  We  have  also  been  exploring  quasi-analytical  models  of  the  ablation  process 
for  fragile  char  systems  to  better  understand  which  material  properties  are  most  important  to  the 
overall  ablation  process.  T  he  issues  associated  with  each  of  these  areas  will  be  discussed  in  the 
following  sections. 

Population  Balance  Modeling  of  Polymeric  Systems 

As  previously  noted  it  is  possible  to  formulate  a  population  balance  based  approximation  for  the 
polymer  degradation  process  as  an  intermediate  modeling  formalism  between  molecular 
dynamics  based  methods  and  continuum  based  models.  The  intermediate  approximation  will  be 
used  to  better  understand  kinetics  of  degradation  and  will  facilitate  development  of  more 
rigorously  formulated  continuum  kinetics.  Madras  and  McCoy  ,  Kruis  et  al.  and  Staggs 
describe  population  balance  equations  for  polymer  scission  and  aggregation  (recombination). 
One  imagines  that  a  polymer  can  be  characterized  as  a  population  of  individual  particles  of 
molecular  weight  v.  The  polymer's  molecular  weight  density  function  is  specified  as  P(v), 
where  v  is  a  continuous  variable  representing  the  polymer  molecular  weight.  In  general,  the 
polymer  elements  can  undergo  scission  and  recombination.  Madras  and  McCoy  discuss  a 
general  mass  balance  statement  for  a  random  scission  that  can  be  written  as: 

P{v)  K*  >P(v')  +  P(v-v')  (1) 

Note  that  the  subscript  of  the  rate  coefficient  Kr  represents  a  random  scission.  Madras  and 
McCoy  27  also  consider  the  possibility  of  a  chain  end  scission  with  rate  coefficient  Ki.  Assuming 
first  order  degradation,  one  can  write  a  general  dynamic  equation  for  the  polymeric  system  as: 
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00 

dP(v,t)  /  dt  =  -Kr(v)P(v)  +  2  j  K  R(v')P(v')Q(v,v')dv' 

V 

00 

-  Kj(v)P(v,t)  +  2  jAr/(v')P(v')Q(v,v')flfv' 

V  (2) 

V 

+  1/  2  \/3(v',v  -  v')P(v')P(v  -  v')dv' 

0 

00 

-  j  /3(v,v')P(v)P(v  ')dv' 
o 

In  the  above  equation  Q  is  the  transition  probability  indicating  how  a  polymer  of  molecular 
weight  v  fragments  into  two  fragments  of  v  and  v-v’.  The  recombination  or  aggregation  kernel 
is  (3 .  The  recombination  kernel  indicates  the  frequency  of  collision  and  aggregation  of  v’  and  v- 
v’  in  forming  a  v  molecular  weight. 

Various  approximations  are  available  in  the  literature  for  the  respective  kernels.  Madras  and 
McCoy'  provide  an  expression  for  the  scission  kernel. 

Q(v,  v1)  =  vm  (v  *-  v)m  T(2  m  +  2)  /  [r(w  + 1)2  (v  ’)2m+1  ]  (3) 

They  note  that  m=0  corresponds  to  random  scission  and  m— >oo  corresponds  to  midpoint  scission. 
For  chain  end  scission,  the  kernel  becomes 

Q(v-v„v')  =  £[v-(v'-v,)]  (4) 


The  non-linear  integral-differential  equation  governing  the  growth  of  the  polymer  is  quite  hard  to 
solve.  The  same  equation  occurs  in  aerosol  systems  and  is  known  there  as  the  General  Dynamic 
Equation  (GDE).  Only  for  very  simple  classes  of  problems  do  analytical  solutions  exist 
(Williams  and  Loyalka30).  Various  numerical  solution  techniques  have  been  put  forward  to  solve 
these  types  of  problems.  We  have  studied  these  problems  in  our  research  (Upadhyay  and 
Ezekoye,31  32  Upadhyay33).  An  obvious  technique  for  solution  of  the  above  equations  is  to 
discretize  the  size  distribution  into  differential  sections  (Gelbard  and  Seinfeld34).  The  so-called 
sectional  thus  requires  sectioning  the  molecular  weight  size  distribution  function  into  a  finite 
number  of  size  bins.  Equations  are  solved  for  the  fragment  number  concentration  in  each  size 
bin.  All  fragments  contained  in  a  size  bin  are  considered  to  have  the  same  size.  The  interactions 
between  the  size  distributions  are  modelled  as  interactions  between  fragments  consisting  of 
discrete  sizes.  There  are  computational  challenges  in  maintaining  the  mass  conservation 
properties  of  the  distribution  when  using  sectional  methods.  Although  recent  advances  have 
eliminated  most  of  the  serious  numerical  errors  associated  with  the  sectional  method  it  is  still 
computationally  expensive  when  dealing  with  problems  involving  spatial  variation  of  the  size 
distribution  function.  Sectional  methods  are  certainly  not  the  best  when  used  in  conjunction  with 
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continuum  computational  transport  models  that  solve  for  the  transport  of  heat  and  chemical 
species. 

The  major  problem  with  the  sectional  technique  is  that  the  solution  attempts  to  fully  reconstruct 
the  size  distribution  function.  There  is  considerable  computational  overhead  associated  with 
such  a  calculation.  To  calculate  most  of  the  physical  properties  of  the  polymer  system  does  not 
require  a  full  knowledge  of  the  size  distribution.  This  is  especially  true  of  properties  most 
relevant  to  engineering  applications  such  as  total  volume  loading,  average  size,  polydispersity 
index  etc.  Some  properties  useful  for  engineering  applications  can  be  defined  in  terms  of  the 
moments  of  the  distribution  function  and  are  shown  in  Table  3.  Conceptually,  moment  methods 
are  based  upon  evolution  equations  for  the  moments  (weighted  averages  of  the  distribution).  The 
equations  governing  the  evolution  of  the  moments  are  generally  simpler  to  solve  than  the 
equations  governing  the  complete  evolution  of  the  size  distribution.  Further,  it  is  possible  to 
incorporate  the  moment  evolution  equations  into  computational  frameworks  which  require 
particle  dynamics  coupled  to  other  transport  dynamics. 

A  moment  is  a  weighted  integral  average  over  the  size  distribution  function.  Populations 
described  by  more  than  one  variable  give  rise  to  mixed  moments.  For  example,  if  particle 
molecular  weight  is  chosen  as  a  descriptor  of  the  particle  then  the  various  (integer)  moments  are 
defined  as 

00 

Mm(x,t)=  jvmP(v,x,t)dv,(m  =  0,1, 2, 3 . )  (5) 

0 

Transforming  the  general  dynamic  equation  into  a  moment  equation  requires  multiplying  the 
equation  by  a  power  law  of  the  molecular  weight  and  then  integrating  the  subsequent  equation 
over  all  molecular  weights.  For  the  simpler  case  of  no  recombination  and  the  forms  of  the 
kernels  specified  earlier,  Madras  and  McCoy  provide  the  form  of  the  moment  equation  as: 


dMj  I  dt-‘ 


-kromj+a 


'  2ZJmKR0M 


•  KI0MJ+r 


'j' 


d= o 


vw/ 


In  the  summation  we  use  the  standard  notation  for  a  binomial  expansion 


\u / 


(6) 


Note  in  the  above 


equation  that  there  were  terms  derived  from  integrals  over  the  unknown  size  distribution  which 
may  not  be  explicitly  expressed  in  terms  of  the  moments  to  be  computed  (Upadhyay33).  To  be 
able  to  solve  the  equation,  the  terms  involving  the  size  distribution  needs  to  be  approximated  in 
terms  of  the  moments  alone.  Hulburt  and  Katz*'  introduced  the  notion  of  a  closure  problem  for 
moment  based  simulation.  To  get  a  closed  set  of  moment  equations,  these  integrals  involving  the 
unknown  size  distribution  need  to  be  approximated  using  a  finite  number  of  moments.  In  cases 
where  this  cannot  be  done  easily,  a  large  number  of  moments  may  be  needed  requiring  the 
solution  of  a  large  number  of  moment  equations.  This  is  not  always  feasible  as  the  higher 
moments  give  more  weight  to  the  tail  of  the  distribution,  which  in  most  cases  is  not  very  well 
known  and  is  also  not  very  useful.  One  way  out  of  this  problem  is  to  assume  a  standard  form  for 
the  size  distribution  function.  Usually  the  lognormal  and  gamma  functions  have  been  used.  As 
the  dynamics  unfold,  the  form  of  the  size  distribution  is  assumed  to  remain  unchanged.  The 
changes  in  the  size  distribution  are  assumed  to  occur  only  due  to  the  changes  in  the  adjustable 
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parameters.  The  equations  governing  the  change  of  these  parameters  are  obtained  by  substituting 
the  assumed  distribution  into  the  GDE  and  taking  the  necessary  moments.  The  application  of 
this  technique  along  with  some  distributions  that  are  physically  realistic  and  mathematically 
simple  is  given  in  Williams  and  Loyalka.30  Since  the  lognormal  distribution  is  a  fairly  good 
representation  for  the  aerosol  volume  distribution  in  most  reactors,  it  has  been  used  most  widely 
in  aerosol  reactor  analyses. 

Since  in  the  moment  methods  we  are  not  concerned  with  the  size  distribution  function,  a  simpler 
and  more  general  method  would  be  to  treat  the  size  distribution  as  an  unknown  weight  and  apply 
Gaussian  integration.  This  is  the  basis  of  the  quadrature  method  of  moments  (QMOM)  which  we 
will  discuss  in  more  detail. 

Table  3.  Physical  Relevance  of  Moments  and  Moments  Ratio. 


Moment 

Physical  relevance 

M„ 

Total  number  concentration 

M, 

Weight  fraction 

M,/M0 

Number  averaged  molecular  weight 

m2/m, 

Weight  averaged  molecular  weight 

The  quadrature  method  of  moments  (QMOM)  is  a  method  for  evaluating  integrals  that  appear  in 
the  General  Dynamic  Equation.  The  unclosed  terms  are  closed  by  using  Gaussian  quadrature. 
The  task  is  to  evaluate  these  integrals  as  accurately  as  possible  when  the  size  distribution  P(v,t) 
is  unknown  but  a  certain  number  of  its  moments  are  known.  We  have  explored  solutions  of  the 
GDE  using  QMOM  and  are  in  a  position  to  solve  for  degradation  process  that  has  been  outlined 
above. 

Quasi-analytical  Modeling  of  Ablation  Parameters  (Property  Sensitivity) 

We  have  also  been  investigating  quasi-analytical  solutions  of  the  ablation  process.  This  is 
instructive  because  it  allows  us  to  perform  sensitivity  analysis  to  better  understand  the  how 
materials  might  be  selected  in  a  general  ablation  process.  Because  ablation  is  a  complex 
phenomenon  involving  multiple  modes  of  insulation  and  surface  recession,  it  is  difficult  to 
design  an  optimal  ablative  material.  Ablative  materials  protect  the  back-wall  structure  through 
conductive  resistance,  endothermic  chemical  reactions,  as  well  as  more  complex  mechanisms 
such  as  pyrolysis-gas  “blowing”.  Also,  the  corresponding  surface  recession  of  the  ablative 
material  reduces  the  ability  of  the  material  to  protect  the  back  surface.  Many  of  these 
mechanisms  of  insulation  and  surface  recession  are  coupled.  Therefore,  improving  the  capacity 
of  an  ablative  material  to  perform  one  of  these  mechanisms  might  hinder  the  capacity  of  that 
material  to  perform  other  ablation  mechanisms. 

In  order  to  better  understand  ablation  phenomena,  it  is  beneficial  to  know  which  ablation 
mechanisms  are  most  influential  in  reducing  the  back-wall  temperature.  The  goal  of  this  study  is 
to  produce  some  rules  of  thumb  that  are  applicable  for  specific  types  of  ablative  materials  (e.g. 
simple,  non-charring,  melting,  etc.)  in  specific  operating  conditions  (e.g.  input  heat  flux,  particle- 
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laden  flows,  etc.).  In  the  following,  we  used  an  analytical  model  of  a  non-charring  ablative 
material  developed  by  Braga  et  al.  6  This  model  will  henceforth  be  referred  to  as  the  Braga 
model.  After  the  model  is  presented  and  our  code  is  validated,  we  perform  a  sensitivity  analysis 
to  study  the  relative  impact  of  four  thermophysical  properties  on  the  time  to  failure  of  the 
protected  material.  It  is  found  that  material  density  is  the  most  important  parameter  for 
increasing  the  failure  time  for  the  type  of  material  and  operating  conditions  studied.  Finally,  we 
conclude  with  some  ideas  for  extending  this  type  of  study  of  ablation  phenomena. 

Approximate  Analytical  Model  of  One-dimensional  Ablation 

For  our  sensitivity  analysis,  we  will  use  the  closed  form,  approximate,  analytical  solution  for 
one-dimensional  ablation  developed  by  Braga  et  al.  The  model  is  developed  with  the  aid  of  the 
integral  method  (Goodman37),  and  assumes  a  simple  ablative  material — that  is,  a  material  that 
does  not  form  multiple  phases  in  the  ablation  process.  The  Braga  model  is  developed  for  a  semi- 
infinite  solid.  For  our  analysis,  this  means  that  we  assume  that  the  back-wall  material  has  the 
same  material  properties  as  the  ablative  material.  Braga  et  al.  have  developed  other  models  that 
assume  that  the  back  wall  is  perfectly  insulated  (Braga38).  We  also  make  the  following 
assumptions:  the  material  properties  are  constant  with  temperature  and  time,  the  input  heat  flux 
is  constant,  and  the  material  is  immediately  removed  from  the  surface  after  reaching  the  ablation 
temperature.  A  schematic  of  the  model  is  presented  in  Figure  4,  below. 


Figure  4.  Schematic  of  non-charring  ablative  material. 


At  this  point,  two  features  of  this  schematic  should  be  emphasized.  First,  the  spatial  domain,  x , 
is  divided  into  three  regions:  removed  material,  heated  material,  and  non-heated  material.  The 
model  described  below  uses  different  modeling  equations  to  describe  the  temperature 


17 


•  • 

Experimental  and  Numerical  Characterization  of  Polymer  Nanocomposites  for  Solid  Rocket  Motor 
Insulation  Materials  Annual  Report  (9-30-06) 


distribution  in  these  three  regions.  The  second  feature  that  should  be  emphasized  is  that 
temperature  is  constrained  between  the  ablation  temperature,  TA,  and  the  initial  temperature  of 
the  material,  T0  .  The  primary  purpose  of  the  model  is  to  describe  how  the  temperature 

distribution  of  the  heated  material  region  varies  with  respect  to  time.  Given  this  temperature 
distribution,  the  boundaries  of  the  three  regions  can  be  determined  with  respect  to  time. 


The  fundamental  equation  describing  one-dimensional,  transient  heat  transfer  is  given  by 


pcP 


dT(x,t ) 
dt 


d  (hdT(x, t)} 
/c 

5xl^  dx  J 


(7) 


The  integral  method  defines  some  penetration  depth  (usually  heat  penetration  depth),  8 ,  over 
which  the  above  equation  can  be  integrated.  Braga  et  al.  define  both  an  ablation  depth,  dA(t), 
and  a  heat  penetration  depth,  5P  (t) .  The  domain  of  integration  is  between  these  two  depths.  If  a 
temperature  distribution  is  assumed,  both  T(x,t )  and  S(t)  can  be  determined  by  a  closed-form 
approximate  solution.  The  Braga  model  specifies  temperature  to  follow 


f 


T(x,t)  =  A(t) 


Sp  (Q  -  x 
KSp(t)-SA(t) 


\» 


+  B(t) 


(8) 


where  Aft)  and  Bft)  are  parameters  used  to  adjust  the  temperature  distribution  to  a  given  set  of 
boundary  conditions.  The  parameter  n  can  be  adjusted  in  order  to  provide  a  more  precise 
temperature  distribution  solution. 

The  Braga  model  defines  two  distinct  time  periods.  In  the  pre-ablation  period ,  the  only  physical 
phenomenon  taking  place  is  heat  conduction  through  the  solid.  Ablation  has  not  begun  because 
at  no  point  is  T(x,t )  >  TA .  Therefore,  the  ablation  depth  is  constant  with  time.  In  the  ablation 
period,  surface  recession  and  heat  absorption  due  to  the  heat  of  ablation  begin.  For  both  of  these 
periods  a  different  set  of  boundary  conditions  are  applied  resulting  in  different  equations  for 
T(x,t),  SA(t),  and  Sp(t) . 

For  the  pre-ablation  period,  the  ablation  depth  is  constant.  That  is, 

dSA(f) 

dt 


Although  the  Braga  model  allows  for  the  case  where  there  is  some  initial  surface  recession 
before  ablation  begins,  we  will  assume  that  the  initial  ablation  is  zero.  The  temperature 
distribution  solution  is  constrained  by  the  following  boundary  conditions: 


dTjxJ)  „ 
dx  q  ’ 


at  x  =  SA  ( t ) 


18 


•  • 

Experimental  and  Numerical  Characterization  of  Polymer  Nanocomposites  for  Solid  Rocket  Motor 
Insulation  Materials  Annual  Report  (9-30-06) 


t  0. 

dx 

T(x,t)  =  T0, 
T(x,t)  =  T0, 


atx  =  SP  (t) 

atx  =  SP(t ) 
&tt  =  tn 


Applying  the  integral  method  to  equation  7  with  the  temperature  distribution  in  equation  8,  and 
assuming  the  preceding  boundary  conditions,  Braga  et  al.  derive  the  following  equations  for  the 
ablation  depth  and  the  temperature  distribution: 


SP(0  =  . 

T(x,t)  = 


I  kn(n  + 1)/^ 


pcF 

q"SP{t)(  5Pjt)-x 

Spit ) 


kn 


Y 


+T0, 


T 

i0’ 


x  <  5P{t) 
x  >  SP(t) 


(9) 

(10) 


Note  that  equations  9  and  10  are  not  in  the  same  form  as  they  are  presented  by  Braga  et  al. 
Specifically,  we  made  several  simplifications  because  we  are  assuming  a  constant  input  heat  flux 
and  no  initial  surface  recession. 


Once  ablation  begins,  a  new  set  of  modeling  equations  becomes  relevant.  The  boundary 
conditions  for  the  ablation  period  are 


T(x,t)  =  TA,  x  =  SA{t) 


dx 

dx 

T(x,t)  =  T0, 


dt 


at  x  =  SA  ( t ) 


at  x  =  Sp  ( t ) 
at  x  =  Sp  ( t ) 


The  last  two  boundary  conditions  are  the  same  as  in  the  pre-ablation  period.  The  second 
boundary  condition  above  now  specifies  that  the  heat  into  the  solid  is  reduced  by  a  heat  of 
ablation  term  that  is  directly  proportional  to  the  surface  recession  rate.  Proceeding  as  before, 
Braga  et  al.  derived  equations  for  heat  penetration  depth,  ablation  depth,  and  temperature 
distribution 
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SA(J)  =  SA 
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(13) 


In  the  preceding  equations,  v  =  A,/(cp(Ta  -  T0 ))  is  the  inverse  Stefan  number.  Again,  equations 
1 1  to  13  are  slightly  simplified  versions  of  the  equations  presented  by  Braga  et  al. 

The  model  described  above  is  sufficient  to  determine  T(x,t )  given  the  four  thermophysical 
properties  ( p,  k,  cp ,  and  X )  and  the  initial  and  ablation  temperatures  ( T0  and  TA,  respectively). 

In  order  to  perform  a  sensitivity  analysis  of  the  four  thermophysical  properties,  we  implemented 
this  model  in-house.  In  order  to  validate  our  code,  we  simulated  the  ablation  of  Teflon  at  a 
constant  heat  flux  of  250  Btul  frs  .  The  resulting  temperature  distribution  and  surface 
recession  rate  were  compared  to  the  results  presented  by  Braga  et  al.  Braga  and  coauthors 
validated  their  model  by  comparing  their  results  with  those  of  Blackwell'  )  Temperature  profiles 
at  times  of  0.2  s,  1  s,  2  s,  3s,  and  4  s  are  presented  in  Figure  5.  The  profiles  on  the  left  are  those 
computed  by  Braga  et  al.  The  points  on  the  left  represent  the  results  computed  by  Blackwell. 
The  profiles  on  the  right  were  obtained  using  our  in-house  code.  It  is  evident  that  our  code 
produces  results  for  the  ablation  of  Teflon  in  close  agreement  with  those  of  Braga  et  al.  and 
Blackwell. 


Pos  ton  (ft) 


Figure  5.  Comparison  of  temperature  profiles  obtained  in-house  (right)  with  results  of 
Braga  et  al.  and  Blackwell  (left). _ 
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As  a  demonstration  of  how  surface  recession  proceeds  with  time,  we  present  plots  of  SA  ( t ) , 
8p{t),  and  u(t)  =  SP(t)~ SA(t)  in  Figure  6.  It  is  seen  that  a  steady  surface  recession  rate  is 
reached  after  approximately  1  s.  After  this  point,  the  ablation  depth  and  the  heat  penetration 
depth  increase  at  the  same  rate.  The  plots  in  Figure  are  also  in  close  agreement  with  Figure  4  in 
Braga.'1' 

In  this  section,  we  have  introduced,  implemented,  and  verified  a  model  for  one-dimensional, 
non-charring  ablation.  We  will  use  this  model  in  the  next  section  to  perform  a  sensitivity 
analysis  of  the  failure  time  of  an  ablative  material  to  four  of  its  thermophysical  properties. 


Figure  6.  Penetration  depths  versus  time.  Matches  Fig.  4  in  Braga. 


Sensitivity.  Analysis 

The  purpose  of  an  ablative  material  is  to  protect  some  structure  from  severe  thermal  attack.  The 
performance  of  an  ablative  can  be  evaluated  in  terms  of  the  time  it  takes  the  back  wall  of  the 
ablative  material  to  reach  some  critical  temperature.  Specifically,  it  is  known  that  the  adhesive 
used  to  bond  the  ablative  to  the  protected  structure  fails  at  approximately  500  F  (960  R). 
Therefore,  we  define  failure  time,  tF,  as  the  time  it  takes  for  the  back  wall  temperature  to  reach 
960  R  after  being  exposed  to  a  heat  flux  at  its  surface. 

In  order  to  study  the  relative  importance  of  the  different  material  thermophysical  properties  to 
tF ,  it  is  useful  to  perform  a  simple  sensitivity  analysis.  To  do  this,  we  first  specified  a  set  of 
nominal  thermophysical  properties  for  Teflon.  These  values  were  obtained  from  Braga.  ’1'  The 
nominal  thermophysical  properties  used  are  p -\20lbm  /  ft2  ,  k  =  3.6xl0~5  Btu!  ft  ■  s  ■  R  , 
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cP  =  0.3  Btul lbm-  R  ,  and  A  -\000  Btu / lbm  .  Running  the  Braga  model  with  these  nominal 
properties  results  in  a  nominal  failure  time  of  tF  =  247.4  s . 


The  sensitivity  of  failure  time  to  a  property  was  determined  by  “perturbing”  the  nominal  value  of 
that  parameter  by  +/-  10%.  For  parameter  x,  we  define  the  positive  “perturbation”  of  x  to  be 
Ax  =  x+  —x  .  For  a  10%  positive  perturbation,  Ax/x  =  0.1 .  Running  the  simulation  with  the 
perturbed  value  results  in  a  perturbed  failure  time,  tF(x+) .  Consequently,  we  can  determine  a 
failure  time  perturbation,  A tF  =  tF(x+)-7F  .  The  sensitivity  of  failure  time  to  a  particular 
parameter  is  then  defined  in  terms  of  the  ratios  of  perturbations.  Specifically,  the  positive 
sensitivity  for  parameter  x  is  found  by 


S 


+ 

x 


A tj_ 

Ax 


(14) 


The  sensitivity  is  normalized  by  the  ratio  of  the  nominal  values  in  order  to  make  a  sensitivity  of 
one  correspond  to  the  case  where  a  10%  change  in  the  parameter’s  value  results  in  a  10%  change 
in  failure  time.  The  negative  sensitivity  for  parameter  x  is  found  by 


At: 


■Ax 


(15) 


Ideally,  positive  and  negative  sensitivities  should  be  approximately  equal.  If  it  is  found  that  this 
is  not  the  case,  it  might  be  necessary  to  perform  a  more  fine-grained  analysis  with  perturbations 
on  the  order  of  +1-5%  or  +/-1%. 

Positive  and  negative  sensitivities  were  computed  for  all  four  thermophysical  properties.  The 
results  of  this  sensitivity  analysis  are  presented  in  Figure  7.  It  is  found  that  failure  time  is  most 
sensitive  to  the  density  of  the  ablative  material.  The  sensitivities  for  p  were  both  approximately 
one.  This  means  that  a  10%  increase  in  density  will  result  in  a  10%  increase  in  failure  time. 
Failure  time  was  least  sensitive  to  thermal  conductivity.  In  fact  a  10%  increase  in  k  resulted  in 
an  imperceptible  change  in  failure  time  for  the  precision  (0.1  s )  of  our  simulation. 

Note  that  the  information  presented  below  gives  no  indication  of  the  direction  of  the  change  in 
failure  time  with  a  given  parameter.  With  the  exception  of  k,  however,  it  was  found  that  failure 
time  was  increasing  with  all  parameters. 
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□  Sensitivity  (+10%) 

□  Sensitivty  (-10%) 


0  0.1  0.2  0.3  0.4  0.5  06  07  0.8  0.9  1  1  1 


Sensitivity  of  Parameter 

Figure  7.  Positive  and  negative  sensitivities  of  failure  time  to  relevant  physical 
parameters. _ 

The  results  of  this  analysis  show,  at  least  for  non-charring  ablatives,  that  conductive  resistance  is 
not  a  relatively  important  mechanism  for  insulating  the  protected  structure.  Rather,  the 
endothermic  reactions  grouped  under  the  heat  of  ablation  contribute  the  most  to  the  protection  of 
the  back  wall  structure.  Increasing  density  means  that  there  is  more  material  to  burn,  and 
increasing  the  heat  of  ablation  means  that  more  heat  is  absorbed  per  unit  mass  of  burned 
material.  While  it  appears  that  specific  heat,  cp ,  is  of  intermediate  importance,  the  primary 
focus  of  the  ablative  material  designer  should  be  increasing  the  heat  of  ablation  capacity  of  the 
material. 

It  should  be  noted  that  these  results  are  preliminary,  and  other  mechanisms  of  insulation  might  be 
found  to  be  more  important  for  different  types  of  ablative  materials.  In  particular,  it  is  desirable 
to  know  how  these  results  will  change  for  a  charring  ablative  material.  Additionally,  the  model 
used  was  chosen  for  its  simplicity.  It  is  not  at  all  clear  how  wide  its  range  of  applicability  is. 
Although  the  Braga  model  provides  accurate  results  for  Teflon,  it  is  not  known  how  well  it 
models  different  non-charring  ablatives. 

5.4  Publications 

The  following  papers  are  presented  or  submitted  for  publications  that  are  sponsored  by  this 
research  program: 

1.  Koo,  J.  H.,  D.  W.  Ho,  and  O.  A.  Ezekoye,  “A  Review  of  Numerical  and  Experimental 
Characterization  of  Thermal  Protection  Materials  -  Part  I.  Numerical  Modeling”,  Paper 
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AIAA-2006-4936,  42nd  AIAA/ASME/SAE/ASEE  Joint  Propulsion  Conference,  Sacramento, 
CA,  July  9-12,  2006. 

2.  Ho,  D.  W.,  M.  Bruns,  J.  H.  Koo,  and  O.  A.  Ezekoye,  “A  Review  of  Numerical  and 
Experimental  Characterization  of  Thermal  Protection  Materials  -  Part  II.  Material 
Properties  Characterization,”  submitted  for  presentation  at  the  48th  A1AA/ASME/ASCE/AHS 
Structures,  Structural  Dynamics,  and  Materials  Conference,  Waikiki,  HI,  Apr.  23-26, 
2007. 

3.  Ho,  D.  W.,  M.  Bruns,  J,  H.  Koo,  and  O.  A.  Ezekoye,  “A  Review  of  Numerical  and 
Experimental  Characterization  of  Thermal  Protection  Materials  -  Part  III.  Experimental 
Testing,”  submitted  for  presentation  at  the  43rd  AIAA/ASME/SAE/ASEE  Joint  Propulsion 
Conference,  Cincinnati,  OH,  July  8-1 1,  2007. 

4.  Nguyen,  K.,  J.  H.  Koo,  D.  W.  Ho,  M.  Bruns,  and  O.  A.  Ezekoye,  “Experimental 
Characterization  of  Thermoplastic  Elastomer  Nanocomposites  under  Simulated  Fire  and 
Solid  Rocket  Motor  Environments,”  submitted  for  presentation  at  the  43rd 
AIAA/ASME/SAE/ASEE  Joint  Propulsion  Conference,  Cincinnati,  OH,  July  8-11,  2007. 

6.0  Summary  and  Future  Plan 

6.1  Experimental 

For  this  period  of  performance,  we  reviewed  numerical  modeling  literature  for  thermal 
protection  materials,  processed  selective  TPUN  formulations  for  material  characterization,  and 
calibrated  a  radiant  heat  panel  apparatus  to  conduct  low  heat  flux  experiments  using  the  TPUNs. 
A  user  proposal  entitled,  “Thermophysical  Properties  Characterization  of  Thermoplastic 
Elastomer  Nanocomposites”  was  approved  by  Oak  Ridge  National  Laboratory/High 
Temperature  Materials  Laboratory  (ORNL/HTML).  This  allows  us  to  use  their  unique  material 
characterization  equipment  and  expertise  to  train  our  graduate  students  to  conduct 
thermophysical  property  measurements  of  our  TPUNs. 

For  the  next  period  of  performance  the  following  tasks  are  planned: 

•  Characterize  selective  TPUNs  using  TGA,  DSC,  TEM,  and  SEM  analyses  at  UT 

•  Prepare  virgin  and  charred  TPUN  specimens  for  specific  heat,  thermal  diffusivity,  and 
thermal  expansion  measurements  at  ONRL/HTML 

•  Conduct  characterization  experiments  using  radiant  heat  panel  (UT)  and  simulated  solid 
rocket  motor  (TSU) 

•  Design  and  fabricate  TPUN  multi-function  sensor  (NANMAC)  for  LHMEL  testing 

•  Conduct  LHMEL  testing  (AFRL/WPAFB) 

6.2  Numerical 

For  this  period  of  performance,  we  summarized  a  simple  model  for  non-charring  ablation,  and 
then  used  that  model  to  perform  a  sensitivity  analysis  of  failure  time  to  four  thermophysical 
properties.  We  found  that  failure  time  is  most  sensitive  to  the  density  and  the  heat  of  ablation  of 
the  material.  It  was  also  found  that  failure  time  is  relatively  insensitive  to  the  thermal 
conductivity  of  the  material.  These  results  indicate  that  improving  the  heat  of  ablation  capacity 
of  an  ablative  material  should  be  the  focus  of  future  research. 
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Because  these  results  are  preliminary,  it  is  necessary  to  confirm  these  trends  using  more  complex 
models — particularly  models  that  take  into  account  the  charring  of  the  material.  For  charring 
ablatives,  it  might  be  that  the  thermal  conductivity  of  the  char  layer  is  considerably  more 
important  that  the  heat  of  ablation. 
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